Change of decoherence scenario and appearance of localization due to reservoir 

anharmonicity 
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Although coupling to a super-Ohmic bosonic reservoir leads only to partial dephasing on short 
time scales, exponential decay of coherence appears in the Markovian limit (for long times) if 
anharmonicity of the reservoir is taken into account. This effect not only qualitatively changes the 
decoherence scenario but also leads to localization processes in which superpositions of spatially 
separated states dephase with a rate that depends on the distance between the localized states. As 
an example of the latter process, we study the decay of coherence of an electron state delocalized 
over two semiconductor quantum dots due to anharmonicity of phonon modes. 



Decoherence of open quantum systems has become one 
of the central issues of the quantum theory. On one hand, 
erasure of phase information with respect to a certain ba- 
sis of states may bring classical behavior out of quantum 
evolution [3- 13- l3| - On the other hand, the ability to con- 
trol the states of quantum systems may open the way to 
novel applications, like quantum information processing 
0, 0, provided that quantum coherence is maintained 
over a sufficiently long time. Therefore, understanding 
and reducing decoherence is also of practical importance. 

Of special interest is a class of models that allow for 
erasure of phase information without transitions between 
the selected basis states. Such a pure dephasing process 
is an essential ingredient of a measurement, with the ba- 
sis states selected by the coupling between the measure- 
ment device and its environment (the pointer basis) jfj 
and determining the physical meaning of the measure- 
ment |3( • A simple model leading to this kind of behav- 
ior is composed of a two-level system and a bosonic bath 
with the coupling between these two subsystems linear in 
the bosonic operators and commuting with the Hamilto- 
nian of the system (independent boson model |?J). Pure 
dephasing effects are also relevant for the short time dy- 
namics of confined carrier states in semiconductor quan- 
tum dots (QDs) H, @ • A class of systems described by 
such a model (including confined carriers in QDs inter- 
acting with phonons) shows only partial dephasing with 
a finite asymptotic level of coherence 0,^3- This feature 
of real systems is essential for the possibility of explain- 
ing the classical nature of measurement results in terms 
of dephasing ( "einselection" ) Q , which is based on the 
expectation that dephasing processes are asymptotically 
exponential, as observed in some formal models |ll| . 

The purpose of this paper is to show that the important 
property of complete or partial asymptotic dephasing de- 
pends not only on the system-reservoir coupling but also 
on the properties of the reservoir itself. In particular, no 
real reservoir may be strictly harmonic, one obvious rea- 
son being the need to restore the equilibrium, i.e., a suf- 
ficient level of ergodicity. Here, this is accounted for by 
including an anharmonic coupling between the bosonic 
modes. This model closely corresponds to real properties 



of solid-state phonon reservoirs, where also parameters of 
the model may be inferred from experiments. 

The combination of system-reservoir coupling with 
reservoir anharmonicity leads to a dephasing effect that, 
for spatially localized systems, may be understood in 
terms of collisional decoherence: Because of their cou- 
pling to the system, reservoir modes undergo a shift of 
their equilibrium positions which depends on the system 
state, forming a coherent displacement field 0,0]. If 
the modes are coupled by anharmonicity the displace- 
ment field acts as a scattering potential for other reservoir 
modes. Since the displacement is state-dependent, each 
scattering event extracts a certain amount of informa- 
tion on the system state, gradually leading to complete 
dephasing. This is the first essential result of the present 
work. The second result is related to the special situation 
when the system is in a superposition of two states corre- 
sponding to distinct positions in real space. In such case, 
scattering of reservoir modes leads to vanishing of the 
coherence between the two distant states. This process 
turns a genuine quantum-delocalized system state into a 
mixture of two classical-like localized states and is hence 
referred to as localization 0,0,0- It is shown, for a 
specific model of carrier-phonon interaction in semicon- 
ductor QDs, that the rate of localization grows with the 
separation between the two spatial locations. This con- 
firms the intuitive expectation that coherence on large 
distances should be more fragile. 

Let us consider the Hamiltonian H — Hq + Hi , where 
Hq is the two-level independent-boson Hamiltonian 

Ho = -\E^ z +Y,^kb{b k + l -a z Y,n{bk + bl k ) 1 (1) 

k k 

with Ft = F_fc, and Hi describes the third order anhar- 
monic coupling between various phonon modes 

Hi = - Wk 1 ,k 2: k 3 Sk 1 +k 2 +k 3 =QA kl A k2 A k3 , (2) 

k 1 k 2 k 3 

where A k = b k + b_ k and the anharmonic constants 
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are symmetric under permu- 



tation of indices. The polarization (branch) index of the 
boson modes is implicit in k. 
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We define the unitary operator 

w = |oxo|® + |ixi|® w, 



(3) 



where W = exp[(l/2) J2 k 9 k b k - H.c], g fe = F k /(hu k ). 
In terms of the new operators f3 k = WfefcW^ = 
&fc + k&z9ki the Hamiltonian i?o is diagonal, 
H = -(l/2)Ea z + Ysk^kfilflk- Using the exact 
diagonalization by the operator W one can find the 
evolution of the non-diagonal element of the reduced 
density matrix (in the original basis), |poi(i)| = 
\p i(0)\exp[~2jdwcoth(hui/k B T)J(uj)sin 2 (ojt/2)/(huj) 2 } 1 
where J(ui) = J2k \Fk\ 2 S(oj — ui k ). For spectral densities 
sufficiently regular at low frequencies, J(u>) ~ ui n , n > 2 
(super- Ohmic reservoirs) and for gapped reservoirs 
(J (to) — around u> — 0), \poi\ reaches a finite asymp- 
totic value (corresponding to partial dephasing) [8J, [l0| . 
On the other hand, if J(lo) ~ u> (Ohmic reservoirs), poi 
decays exponentially for long times [l7l |. It turns out, 
however, that Ohmic independent boson models show 
infrared divergences that cause problems on the formal 
level ^3 • Below we will see that dephasing generated by 
a super-Ohmic reservoir becomes complete (exponential 
at long times) if the reservoir is anharmonic. 

In terms of the transformed operators (3 the anhar- 
monic Hamiltonian H\ becomes 



Hi — - *^2, w k 1 ,k 2 M^k 1 +k 2 +k 3 =0 (Ak 1 Ak 2 Ak 3 

klk 2 k 3 

-J>o z g kx A k2 A k3 + ^gk 1 gk 2 Ak 3 - crzgk 1 gk 2 9k 3 ) , 

where Ak = /3fe + ft-u- The first term is the anharmonic 
coupling between the new modes, the third describes a 
shift of the oscillator equilibria, the fourth is a shift of 
the energy levels that may be included in the energy E. 

Of interest here is the second term that describes two- 
phonon absorption, emission and scattering. Only the 
latter may lead to energy-conserving processes that do 
not involve real transitions between system states. Us- 
ing commutation relations and symmetries of the anhar- 
monic coefficients one may write the relevant (scattering) 
part of the anharmonic Hamiltonian 



H 



-&z ^ W q -k,-g,k9q-k 



(4) 



k.q 



x [/3 q /3 k - 8 qtk n k + Sg ik (n k + 1/2)] 



The last term may again be included into the energy E 
(it vanishes in the limit of infinite reservoir volume). It 
should be noticed that Eq. 10} is non-diagonal in phonon 
modes which allows for real pho non scattering processes 
(unlike the model of Ref. 18]). A similar scattering 
Hamiltonian may be obtained by including higher exciton 
levels into a purely harmonic model |19| . 

It is known that scattering on a heavy Brownian parti- 
cle leads to decoherence and localization of the quantum 



state of the latter 0, 0, 0] . In order to see that the 
same is true in the present case of scattering of bosonic 
modes on a two-state quantum system and to extract the 
corresponding long-time behavior we write the evolution 
equation with the interaction Hamiltonian . Assuming 
that the time scales of the dephasing process are longer 
than those related to the reservoir memory and initial 
dephasing one may consistently describe the long time 
dynamics in the Markov limit. From the resulting Lind- 
blad equation one finds the solution describing exponen- 
tial pure dephasing at long times with the rate |20| 



1 

T 2 



dt(B^(t)B(0)) 



(5) 



= 4%22 \ W q-k-q,k\ 2 \gq-k\ 2 nk{n q + l)5(u! q - U) k ), 
k.q 

where B(t) = ^ k q w^-^-h^h^'-^* - 
5q, k n k ] are reservoir operators with vanishing equilibrium 
average. The dephasing rate is finite at T > 0, which 
should be contrasted with the harmonic case, where the 
dephasing effect vanishes in the Markov limit. Thus, 
reservoir anharmonicity has qualitatively changed the de- 
coherence properties of the system. 

The rate given by Eq. (0 is consistent with the scat- 
tering picture described in the introduction. A boson 
scatters from the state k to q. The scattering ampli- 
tude is proportional to the magnitude of the displacement 
field, governed by the system-reservoir coupling constants 
g q -k, and to the anharmonic coupling w q — kt —q tk . The 
momentum transfer in such an event cannot exceed the 
inverse size of the displacement field which will be re- 
flected by a cutoff in g q - k (see below). The scattering 
probability depends on the occupations of the initial and 
final states. Finally, since no real transitions between 
system states are allowed, the scattering must be elastic, 
as expressed by the energy conserving Dirac delta. 

In order to see if this effect may be of importance un- 
der realistic conditions, let us now study the specific case 
of a single electron in a pair of semiconductor quantum 
dots, as in the recent coherent manipulation experiment 
[2l| . The decay of coherence between the localized states 
corresponds in this case to a localization process, hence 
Ti may be referred to as the localization time. Since we 
are interested in the dephasing due to anharmonicity, we 
disregard phonon-assisted tunneling between the states 
(cr Xt y coupling) which might appear only if the states 
overlap (such terms obviously lead to exponential deco- 
herence; the anharmonicity effects for such real transition 
processes were studied elsewhere [22,. '23]). The Hamilto- 
nian is therefore 



#QD - 



^ + T,ft(bk + bl k ) 



Y^h^kb k o k , 
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where |0), |1) are the basis states (each localized in one of 
the two dots), e, are the energies of the two states, b\ are 
phonon operators (with respect to the unperturbed equi- 
librium), and ft are coupling constants. This Hamil- 
tonian is transformed to the form of Eq. by the 
canonical transformation (shift) of the phonon modes, 

bk = b k + (fl 0) + fj c 1) )/(2hujk). The effective system- 



reservoir coupling constants are then Fk = fjp — '■ 
This shift modifies also the anharmonic Hamiltonian 
Hi but, since the transformation is independent of sys- 
tem state, no extra coupling will appear. The new lin- 
ear and quadratic terms in Hi may be removed by re- 
diagonalizing the phonon Hamiltonian, which produces 
negligible higher order corrections to the couplings. 

In polar semiconductors, the strongest lattice displace- 
ment (polaron |7|) is related to longitudinal optical (LO) 
phonons, which are also subject to strong anharmonic 
coupling to acoustic phonons 0> 113 • Therefore, the 
present discussion will be restricted to the scattering of 
longitudinal acoustic (LA) phonons on the LO displace- 
ment field. We assume dispersionless LO modes with fre- 
quency O = 54 ps -1 (all values for GaAs 26] ). The elec- 
tron wave functions will be modelled by isotropic Gaus- 
sians of size L and the QDs will be displaced by a distance 
D along z 0| ■ The physical coupling constants are 



Jk 



(04) 
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where e is the electron charge, V the normalization vol- 
ume of the phonon modes, Eq the vacuum dielectric con- 
stant, and e = 70 the lattice part of the relative dielectric 
constant. The Gaussian momentum cutoff reflects the 
momentum conservation and momentum-position uncer- 
tainty for an electron wave packet of size L. 

In an anharmonic process, an LO phonon interacts 
with two LA phonons with linear dispersion ujk = ck up 
to the Debye wave vector kjy = 11 nm -1 , where c = 5150 
m/s is the speed of sound. For this process, the general 
form of the coupling is w q ^k,q,k = (wa/VV)y/qk |27j . 
where q, k pertain to the LA phonons and we neglected 
the dependence on the LO phonon momentum in the nar- 
row range of its relevant values k < 1/L -C ferj. Using 
the measured lifetime To = 9.2 ps of the LO phonon at 
k = 24] one finds, using the Fermi golden rule with 
the anharmonic Hamiltonian J2J), u 
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FIG. 1: Left: the dependence of the localization time on 
temperature for L = 4 nm (inset shows the low-temperature 
sector). Right: the dependence on the distance between the 
states for L — 4 nm (solid lines) and L = 8 nm (dashed). 



Eq. one finds after some algebra J2j 

l_ _ 64 e 2 (k B Tf f°° dx _ x 2 
T 2 ~ 7T 2 T Q h e n 5 e ec J x & 

( XXB 



sin ax 



(8) 
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where xrj = (hck^) / \k^T) , a = y/2D/L, and 



cj>{x) 



duu 



(e u - l) 2 



Substituting the above result along with Eq. JJJ into 



The dependence of the dephasing time Ti on temper- 
ature and distance between the basis states is shown in 
Fig. n Note that at moderate temperatures T 2 is much 
longer than the initial dephasing and reservoir memory 
times (~ 1 ps) so that the results are consistent. For low 
temperatures the localization rate is limited by the num- 
ber of occupied initial states and the resulting allowed 
final states (due to energy conservation). As a result, 
one finds T2 ~ T~ 7 and the effect is extremely weak 
for sub-Kelvin temperatures (see Fig. QJi). However, al- 
ready at T ~ 20 K the dephasing rate is of the order 
of nanoseconds and becomes comparable to typical co- 
herent manipulation times on such structures. At ~ 100 
K the dephasing time drops to several picoseconds even 
for closely spaced dots. This temperature dependence 
is much stronger than in the quantum Brownian motion 
§,[2^. This is not astonishing, since this feature depends 
on the reservoir density of states and the physical nature 
of the coupling so that no universality can be expected 
here. In fact, similar strong temperature dependence has 
been found for localization processes due to scattering of 
light on dielectric balls and on free electrons Q. 

At all temperatures the localization rate is increased 
by over an order of magnitude when the distance D be- 
tween the dots grows from nanometers to micrometers. 
This distance dependence is shown in detail in Fig. \})p. 
It turns out that the dephasing rate grows rather fast 
(~ D 2 ) as long as the wave functions overlap. This de- 
pendence is a general feature in the regime of "ineffec- 
tive single scattering event" 0]- When the states get 
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separated, the increase of the localization rate continues, 
although it becomes only logarithmic. This is related to 
the crossover to the regime of spatially distinct states 
where a single scattering event is sufficient to extract the 
position information. The dependence on the dot size L 
becomes less and less important as the distance grows, 
as shown by the comparison between the dot sizes of 4 
and 8 nm. This should be contrasted with the harmonic 
model [i^ , where the dependence on the separation sat- 
urates while that on the size does not. Therefore, the 
present results cannot be fully explained by merely in- 
voking the known effect of increased carrier-phonon cou- 
pling for separated carriers. The logarithmic asymptotic 
behavior is unexpected on the grounds of a general dis- 
cussion [pol flfij l ; here it results from the long range na- 
ture of carrier-phonon interaction, manifested by the long 
wavelength singularity in Eq. ©. 

The presented results show that anharmonicity of a 
super-Ohmic bosonic reservoir leads to a qualitative 
change in the dynamics of an open system. Exponential 
(Markovian) pure dephasing appears even though only 
partial decoherence was present without anharmonicity. 
For typical coupling properties, the dephasing rate de- 
pends very strongly on temperature. In the case of a 
superposition between two spatially separated states the 
rate of the resulting localization grows also with the spa- 
tial distance between the two states. 

It should be stressed that the dephasing (localization) 
mechanism described here is inherent to physical proper- 
ties of the system and appears universally for any lo- 
calized states embedded in a translationally invariant 
bosonic reservoir with a certain degree of anharmonic- 
ity. As such, it sets material-dependent limits to system 
coherence, independent on any design improvements that 
might eliminate noise sources that dominate dephasing in 
the current experiments |2l| . 

The present result is of importance to a few areas. 
First, it describes an additional dephasing mechanism 
that must be taken into account both in design of de- 
vices relying on quantum coherence and in interpretation 
of experiments. Apart from the localization effect dis- 
cussed here, the anharmonic scattering mechanism will 
contribute, e.g., to the broadening of the zero-phonon 
line in QD spectroscopy pffll Isij . Second, the depen- 
dence of dephasing on the distance in space may affect 
scalability of quantum computing schemes and applica- 
bility of concatenation techniques used in quantum fault- 
tolerant architectures • Third, appearance of dephas- 
ing in models with non-singular, super-Ohmic coupling 
may be of importance for emerging of classicality from 
quantum evolution 
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Derivation of the T 2 time 
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This auxiliary material explains how the long time scale behavior of the system was extracted in 
the Markov limit and gives the technical details of derivation of the final formula [Eq. (8)] for the 
dephasing times in the QD localization example. 
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The system evolution in the Markovian limit is found following the usual approach [? 
(discarding the last term) is written in the interaction picture in the form 



The Hamiltonian H[ 



where B{t) is the reservoir operator defined below Eq. (5). We assume that the initial state of the system is separable 
in terms of the shifted operators 

p{0) = ps(0) ® PBj 

where pB is the thermal equilibrium state of the reservoir. Physically, this corresponds to preparing the system in 
such a way that the polaronic displacement field is formed by adiabatically following the system evolution. This 
assumption corresponds to typical timescales in the system of two QDs used as an example in the paper but, in any 
case, the long-time behavior cannot depend on the way the initial system state was prepared. 

The evolution equation for the density matrix of the total system is written in the usual integral form and the 
reduced density matrix of the two-level system is then obtained by tracing out the reservoir degrees of freedom, 



ps(t) = ~ / dr Tib Hi 



H[ S \p(T) 



where we used the fact that the operator B(t) has a vanishing average at thermal equilibrium so that the 1st order 
term vanishes. 

At this point, a series of standard approximations is introduced, leading to the Born-Markov equation for the 
evolution of the two-level system. First, one assumes that the perturbation of the reservoir state due to the joint 
evolution of the two subsystems is weak. It may then be assumed that at each instant of time the system state is 
separable, with the reservoir at thermal equilibrium, p(s) — ps (s) <g) pe , which leads to the Nakajima-Zwanzig equation 



Ps(t) 



dr Ttb 



H 



H[ s) , Ps (t) ® p B 



dr Re(B(t - r)B(t)) (ct z ps(t)ct z - p s (r)) 



where we use the fact that, for a hermitian operator B(t), one has (B(t)B(t')) = (B(t')B(t))* . Next, one assumes 
that the evolution of the two-level system is very slow and its state changes negligibly during the reservoir memory 
time tm over which the correlation function (B(t)B) is essentially nonzero. This allows us to replace ps{i~) by ps(t). 
Finally, at times t > tm the lower integration limit may be set to — oo. Using the fact that pb is a stationary state, 
the Born-Markov evolution equation for the reduced density matrix of the two-level system may be written as 



where 



Ps(t) = \p [<r z ps(t)(T z 



7 = 4 Re / (B(t)B)dt = 2 



Ps(t)}, 



(B(t)B)dt. 



It is easy to see that the diagonal elements of ps(t) remain constant, while the non-diagonal ones decay with the time 
constant T2 = I/7. The final formula of Eq. (5) follows in a straightforward manner by explicitly calculating the 
correlation function using the definition of the operator B{t), as given below Eq. (5) in the paper. 

This result may be directly used to calculate the T2 time constant for the localization process used as an example 
in the paper. Inserting the coupling constants g q ~k and w q -k.- q ,k into Eq. (5) and changing both summations into 
integrals according to the usual formula (for normalization volume V), J2k ~ > ^ / d 3 k/(2ir) 3 , one gets 



7 T 2 7r 4 r O 5 fe e 



d 3 k / d 3 q 



kq 



(k - qf 



-(k-q) 2 



(k z - q z )D 



8{ck — cq)nk{n q 
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Let us introduce the new variable p = k — q. Next, we express q in spherical coordinates (q, 0, </>) in the reference 
system with the z' axis along p and p itself in spherical coordinates (p, ip) in the original (absolute) reference frame. 
In the new variables the Dirac delta is transformed into 

1 S(p+ 2qsm9) 

S(ck-cq) = -6(q- \q + p\) = 



The result becomes 



c c|sint/| 

32e 2 c 5 r , , r /2 .„ „ r . , r /2 



d-d cos j? 



, , y , - / dqq 2 / dflcostf / dpp 2 I 

-K z T Whe S J q<kD J_ n/2 J\q+p\<kv J-n/2 

q\p + q\ -i(Lp) 2 • 2 P Dsin ^ s (p + 2 i sin0 ) i , i \ 

X ^^ 6 2( ' Sm — 2 ^in^^^ + 1} ' 

where the integration with respect to <j>, <p has been performed. It is now possible to integrate over p and d. Introducing 
the new variables 

hcq V2k B TLx . „ 
x = - — — , t= sin 

k B T he 



we can write the result in the form 



64e 2 c 4 fk B T\ 5 f XD , . e x rV2k D L X / XD df 



sin at 



7r 2 Toil 5 feo£ V ^ c / Jo ( ea: — I) 2 7o i \ at 

where xd = hcku/(k B T) and a = \/2D/L. This formula may be integrated by parts to yield 
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^ ^ / _ sinoaA _ ^ cfe^ / sinaaA 

a; V ax J J x \ ax J \v2fc D i 



Since k^L » 1 the upper limits of both integrals may be extended to oo, which leads to the form of Eq. (8). 



